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Abstract. This paper presents experimental data on single droplet combustion of 
decane in microgravity and compares the results to a numerical model. The primary 
independent experiment variables are the ambient pressure and oxygen mole fraction, 
pressure, droplet size (over a relatively small range) and ignition energy. The droplet 
history (JD 2 history) is non-linear with the burning rate constant increasing throughout 
the test. The average burning rate constant, consistent with classical theory, increased 
with increasing ambient oxygen mole fraction and was nearly independent of pressure, 
initial droplet size and ignition energy. The flame typically increased in size initially, 
and then decreased in size, in response to the shrinking droplet. The flame standoff 
increased linearly for the majority of the droplet lifetime. The flame surrounding the 
droplet extinguished at a finite droplet size at lower ambient pressures and an oxygen 
mole fraction of 0.15. The extinction droplet size increased with decreasing pressure. 

The model is transient and assumes spherical symmetry, constant thermo-physical 
properties (specific heat, thermal conductivity and species Lewis number) and single 
step chemistry. The model includes gas-phase radiative loss and a spherically 
symmetric, transient liquid phase. The model accurately predicts the droplet and 
flame histories of the experiments. Good agreement requires that the ignition in the 
experiment be reasonably approximated in the model and that the model accurately 
predict the pre-ignition vaporization of the droplet. The model does not accurately 
predict the dependence of extinction droplet diameter on pressure, a result of the 
simplified chemistry in the model. The transient flame behavior suggests the potential 
importance of fuel vapor accumulation. The model results, however, show that the 
fractional mass consumption rate of fuel ^ the flame is close to 1.0 for all but the 
lowest ambient oxygen mole fractions. ^ ^ 
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1. Introduction 

Theories and experiments involving single droplet combustion date back to 1953 [1], 
with the first microgravity work appearing in 1956 [2]. The problem of a spherical 
droplet burning in an infinite, quiescent, microgravity environment is a classic problem 
in combustion research with the simplified solution appearing in nearly every textbook 
on combustion (e.g. [3]). This simplified solution predicts that the decrease of the 
droplet diameter squared with time is a constant (d-squared law) ( k ). The solution 
also predicts that the ratio of the flame diameter ( Df ) to the droplet diameter (D) is 
constant and independent of time or droplet diameter. 

Because of the ideal geometry and one dimensional nature of the problem, 
spherically symmetric droplet combustion has been extensively studied. The 
microgravity environment offered by ground-based facilities such as drop towers and 
space-based platforms, is ideal for studying the problem experimentally. Indeed as a 
recent review by Choi and Dryer[4] shows, significant advances in droplet combustion 
have been made by studying the problem experimentally in microgravity and comparing 
the results to one dimensional (spatially) theoretical and numerical treatments of the 
problem. 

The classical analysis makes several simplifying assumptions, including a quasi- 
steady gas-phase, constant thermo-physical properties, unity Lewis number, and 
infinitely fast chemical kinetics. These assumptions provide valuable qualitative 
insight into the physics of droplet combustion but make quantitative agreement with 
experiments difficult. The classical analysis can produce reasonable estimates of k (with 
properly chosen thermo-physical properties). Estimates of Df, however, often differ from 
experimental values by a factor of five or more. Further, the classical analysis cannot 
predict flame extinction. Current theoretical (e.g. [5]) and numerical (e.g. [6]) models 
relax many of these assumptions and show better agreement with experimental data. 
To fully exercise the capabilities of these models and improve them, researchers require 
access to single droplet combustion data over a range of fuels and ambient conditions. 

In this article, we present the results of single droplet combustion experiments and 
compare the results to a numerical model. The fuel for all of the experiments was 
decane. The independent experiment variables were the ambient pressure, oxygen mole 
fraction, initial droplet size, and ignition energy. The model is transient and includes 
single-step, finite-rate chemistry and flame zone radiative loss. 

2. Experiments 

2.1. Experimental Apparatus 

Because complete details of the experimental hardware are available elsewhere[7] only 
a brief description is provided here. A 125 or 230 pm fiber with a small bead (« 2x the 
fiber diameter) supported droplets with an initial diameter between 0.9 mm and 1,8mm. 
A stepper motor driven syringe dispensed the fuel droplets onto the ends of the fibers 
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in normal gravity just before the drop. A small coiled hot-wire, withdrawn immediately 
after ignition, ignited the droplets in microgravity. The ignition duration was typically 
on the order of 1 s. During this time the battery delivered approximately 40 J of energy 
to the hot wire. A simplified thermal analysis of the hot-wire shows that approximately 
10 percent of this (i.e. 4 J) actually heated the gas- phase surrounding the droplet during 
the Is ignition time. The rest of the energy heated the wire or radiated to the chamber 
walls. For reference, the total energy released in burning a 1.5mm decane droplet is 
approximately 60J. 

The atmospheres were blends of air} and nitrogen, mixed by partial pressures. We 
estimate the oxygen mole fraction (X 02 ) to be ±0.002 of the stated value, and the 
pressure ( P ) to be accurate to 1 percent of the indicated value. The uncertainty in 
the atmosphere comes from uncertainty in the chamber pressure, small leaks in the 
combustion chamber, and errors that result from the finite volume of the gas supply 
tubing. 

The data for experiments were video recordings from two orthogonally located, 
black and white CCD video cameras. The first camera provided a magnified backlit 
view of the droplet to obtain the droplet regression history. The second view was of the 
flame and had a small 15 pm diameter SiC fiber in the field-of-view (passing through 
the center of the droplet). Radiant emission from this small fiber helped determine the 
presence (or lack thereof) of the flame in ambient conditions in which the CCD camera 
could not image the dim flame. 

A microcomputer equipped with a black and white frame-grabbing board analyzed 
the video data from the experiments. The droplet diameter in this article is the size 
that results from equating the measured volume or projected area of the droplet to 
that of the equivalent sphere or circle, respectively [8]. In reporting the burning rate of 
the fuel droplet, we present both an average, (k) and an instantaneous ( k ) burning rate 
constant. The k comes from a linear fit of the data between t = O.lt*, and 0.9 4, where 
is the total burn time (igniter withdraw to extinction or burnout) and t is the time (t) 
normalized by the initial droplet diameter squared (D%). The k comes from a mo d ified 
cubic spline fit routine applied to the experimental data. In this method, the droplet 
history ((D/ Do) 2 as a function of i) is divided into equally-sized discrete intervals (At). 
A best-fit third order polynomial is found for each interval, subject to the constraint 
that at each node point there is continuity in (D/Do) 2 , its first (k) and second ( k ) 
derivatives. The k at any time is then the derivative of the best-fit polynomial for that 
particular interval. We found that this provides an excellent fit to the experimental 
data. 

t The air was a precision gas blend of oxygen (0.21 mole fraction) and nitrogen mixed gravimetrically 
by the manufacturer. 
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2.2. Experimental Results 

Figure 1 shows the burning history for a single 1.7 mm decane droplet in a P — 
!90mmHg, X 02 — 0.17 ambient. This test exhibits many features typical to all of 
the tests and is also one of the base cases used in the refinement of the n um erical model 
(described later). The droplet burning history is non-linear. There is a small decrease 
in k initially which is probably due to the presence of the igniter. The igniter is a 
heat source which provides extra energy for vaporization and it is reasonable to expect 
a small transient associated with the withdraw. After this short decrease, k increases 
slowly from 0.5 to 0.6 mm 2 /s until the middle of the test. This increase could be the 
result of a number of transient phenomena in droplet combustion, such as transient 
heating of the liquid phase[9] and/or fuel vapor accumulation [24], The magnitude and 
time-scale of the increase are consistent with both phenomena. The k then increases 
slightly until the end of the burn, where it increases more quickly until extinction or 
burnout. 



Figure 1. Combustion history of a decane droplet initially 1.7 mm in diameter 
burning in a X 02 = 0.17, 190 mmHg ambient. The graph shows a sparse subset 
of the experimental data for clarity. 

The flame behavior in Figure 1 is characteristic of that in most of the nitrogen- 
diluted tests. The flame diameter ( Df ) initially increases in size, plateaus for a 
significant period of time in the middle of the test, and then decreases until the end 
of the test. After reaching the maximum value, the flame responds to the shrinking 
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droplet, although the droplet size decreases 25 percent before the flame responds. 

Finally, Figure 1 shows that the flame standoff ratio (Df/D) continuously increases 
from ignition until extinction. This increase was characteristic of all the nitrogen-diluted 
tests. The rate of increase is nearly constant (linear slope) until late in the burn. Close 
to extinction or burnout, the standoff ratio increases faster. 

Figure 2 shows the effect of ambient oxygen mole fraction at P = 190 mmHg 
and D 0 = 1.6 mm (±0.1 mm) for each test. Both the droplet and flame behavior are 
qualitatively consistent with the classical analysis of single droplet combustion. The 
values of k for the X 02 = 0.15,0.17 and 0.19 tests are 0.52,0.57, and 0.65 mm 2 /s , 
respectively. The relative dependence on ambient oxygen mole fraction is very close to 
that predicted by classical droplet combustion theory [3]. The flame size increases with 
decreasing ambient oxygen mole fraction. All three tests show the same qualitative trend 
of the flame size increases, reaching a maximum, and decreasing until burnout. The 
flame standoff also increases with decreasing ambient oxygen mole fraction, although the 
differences are small between the two smallest ambient oxygen mole fraction ambients. 
The flame standoff increases almost linearly with time for all of the tests. 



Figure 2. Combustion histories for three experiments (Do s=s 1.6mm) burning in 
190 mmHg, and X 02 = 0.15, 0.17 and 0.19 ambients. The graph shows a sparse subset 
of the experimental data for clarity. 


Figure 3 shows the influence of ambient pressure on the burning history of four 
droplets ( D 0 ps 1.7 mm ± 0.1mm) burning in ambients with X 0 2 — 0.15. The droplet 
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regression histories for these tests are nearly identical, with the variation in k (and k ) 
less than 10 percent of the average value of 0.52mm 2 / s. The k starts at approximately 
0.46mm 2 /s and increases throughout the droplet lifetime to approximately 0.65mm 2 /s 
at the end of the test. At pressures below 190 mmHg, the flame extinguishes at a finite 
droplet size. The extinction droplet diameter, D^, increases with decreasing pressure. 
At OOmmHg, the droplet burns for only a short time before flame extinction. 



Figure 3. Combustion histories for four experiments (Do 1.7mm) burning in 
X 02 = 0.15 and P = 760,380,150 and 100 mmHg ambients. The droplet histories 
were nearly identical for all of the tests, so only one is displayed. The graph shows a 
sparse subset of the experimental data for clarity. 


While the droplet behavior was nearly independent of pressure, Figure 3 shows that 
at pressures below 380 mmHg, the flame size (at a given t) increases with decreasing 
pressure, as does Df/D. These increases are consistent with the effects of finite rate 
chemical kinetics that slow the chemical reaction rate down (increasing characteristic 
chemical time) with decreasing pressure. 

Figure 4 shows the effect of initial droplet size for droplets burning in a 100 mmHg, 
0.15 oxygen mole fraction ambient. The flames surrounding droplets burning in this 
ambient pressure and oxygen mole fraction extinguish at a finite droplet size. The 
results show that k did increase with increasing droplet size. The change, however, was 
less than 10 percent over the range of droplet sizes. As one expects, the flame size 
increases with increasing droplet size. The flame standoff dependence on droplet size 
is not as intuitive. The two larger droplets show a consistent, nearly linear increase 
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over the entire droplet burning history. The flame standoff for the smallest droplet, 
however, increases at a slower rate than the two larger droplets, and then reaches an 
almost constant value for the last one third of the burn time. Finally, Fig. 4 shows that 
D ex t (he- extinction droplet diameter) increases with increasing D 0 . 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

t/D a 2 (s/mm 2 ) 


Figure 4. Combustion histories for three experiments with different initial droplet 
diameters burning in Xoi = 0.15 and P = 100 mmHg ambients. The droplet regression 
histories are similar for all of the tests, so only one (Do = 1.54 mm) is displayed 
(k « 0.56 mm 2 /s). The droplet continues to vaporize after flame extinction, thus the 
continued droplet history after the flame data ends. The graph shows a sparse subset 
of the experimental data for clarity. 

Figure 5 shows the effect of ignition energy on the burning history. The variation 
in ignition energy results from changing the time duration of the igniter (r ig ). The two 
ignition times in Figure 5 are 0.7 and 1.2 s and D 0 is 1.7 mm for both tests. Time zero 
in this (and all graphs) corresponds to igniter withdraw. For the longer ignition time, 
a flame appears well before igniter withdraw. The results show that changes in the 
ignition energy do not noticeably affect k, k and minimally affect D ext in these ambient 
conditions. The ignition energy, however, significantly influences flame size and standoff. 
Specifically, this increase in ignition time increases the flame diameter immediately after 
ignition by approximately 25 percent. This difference in flame size persists throughout 
the test with the flame size at extinction much larger for the longer ignition time (higher 
energy) case. 




Single droplet combustion of decane 


8 



Figure 5. Combustion histories for two droplets (Do = 1.7mm) with different ignition 
energies burning in Xoi — 0.15 and P = 90 mmHg ambients. The droplet histories 
were nearly identical for all of the tests, so only one (T ig = 0.7) is displayed. The 
droplet continues to vaporize after flame extinction, thus the continued droplet history 
after the flame data ends. The graph shows a sparse subset of the experimental data 
for clarity. 


3. Numerical Model 
3.1. Model Description 


The numerical model is one-dimensional and transient in both the liquid and gas 
phase. The gas-phase model assumes one-step, second-order overall Arrhenius reaction, 
constant specific heats and thermal conductivity, constant Lewis number for each 
species (although different species can have different, constant Lewis numbers), ideal gas 
behavior, and no buoyant forcefll]. The last assumption allows a simplified treatment 
of the momentum equation[12], which assumes potential flow and that the product (pT) 
is a constant. The igniter is a source term in the energy equation at a particular region 
in the gas-phase. The model includes flame radiative losses from carbon dioxide and 
water vapor (gray gas treatment). The dimensional equations for the gas- phase are: 


l_d_ 2 5£ = Q wp — g r 

r 2 dr r dr PooCpT^ 
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The one-dimensional energy equation for the liquid phase is: 
2 dT . d 2 T l _dT 

a L dt 


r dr ^ dr 2 


( 3 ) 

( 4 ) 

( 5 ) 


The gas-phase radiative loss is q r = Aa{T 4 — T^J. The mean absorption coefficient, 
A, is A = 0.4(P C o2Ap° 2 + T//20^P°)[13]. The multiplication factor of 0.4 reflects the 
non-optically thin nature of the flame and the possible overestimate of the Planck-mean 
absorption data[14]. 

The boundary conditions at the gas/liquid interface are: 


dT dT 

— Ax,— |z, + ^g~Q^\g - Qrs = pu T C 

( 6 ) 

pV F ^ = pu r (l-Y F ) 

( 7 ) 
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pV i ~ = pu r Y i (i = 02, C02, H20 ) 

dr 

( 9 ) 


Equations 6, 7 and 8 represent three equations for the three unknowns at the droplet 
surface, T, Y F and u r . The radiative loss from liquid surface is q rs = ecr(Tj — Tff) where 
e = 0.8. The fair-field boundary conditions are ambient temperature and oxygen mole 
fraction, and no carbon dioxide, fuel or water vapor. The far-field boundary condition 
for the reduced momentum equation is u r = 0 and the boundary condition at the droplet 
surface comes from solving Eq. 6 for u T . 

The fuel is decane and its physical properties are in the nomenclature section. The 
activation energy, E, of the reaction is 30 kcal/mole. The pre-exponential factor, A , 
has a value such that the flame extinguishes when the droplet is at an experimentally 
measured diameter (D ext ), currently 8 x 10 14 cm 3 /g/s). 

^ Pent - 

3.2. Model Solution Procedure 


The governing equations above become a tri-diagonal system of equations after 
discretization. The diffusion terms use a central difference operator while the convective 
terms use an upwind difference. There is an additional convective term to account for 
the fact that the grid moves with time. This is a result of non-dimensionalizing by 
R(t) [9]. The unsteady terms use ' 'a backward difference operator making the entire 
formulation implicit. The non-dimensional radial coordinate is f = tan f, so the domain 
from the droplet surface to the ambient becomes 7t/4 to ir/2. There are 500 variably- 
spaced grid points in the r-direction. Near the droplet surface the smallest cell size is 
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(O.Olr), (O.lf) in the flame zone, and then expands with increasing distance from the 
droplet. The small grid size is necessary near the droplet surface to accurately determine 
the temperature gradient at the droplet surface (which in turn determines the burning 
rate constant and convective field). 

The initial condition is a cold ambient profile with no fuel vapor. The igniter is a 
source term in the energy equation, qj, 4-5 mm from the droplet surface, 1-2 mm in 
width with an energy generation rate between « 1 - 2 cal/cm 3 /s. These values closely 
approximate the energy supplied to the gas phase by the hot-wire in the experiments. 
The source term reduces linearly in magnitude (to zero) as the gas phase temperature in 
the ignition zone increases from 600 K to 1800 K. This decrease simulates the reduced 
heat transfer from the hot wire to the gas phase as the gas temperature increases. 
The results show that the best agreement between the model and experiment (post- 
ignition) occurs when the numerical model closely predicts the experimentally observed 
pre-ignition behavior (droplet size change, ignition delay time). 

3.3. Model Results 

Figure 6 shows a comparison of an experiment and a model prediction for a droplet 
burning in a P = 120 mmHg, X 02 — 0.15 oxygen mole fraction ambient. The D ext are 
0.76 and 0.69 mm for the experiment and model, respectively. The agreement between 
the model and experiment for the droplet burning history is excellent. The predicted 
flame diameter ( D f ) and(1Vs)temporal behavior are consistent with the experiment, 
except near the end of the test. Near extinction, the model predicts that the flame size 
and standoff ( Dj/D ) will decrease until extinction and the experiment shows that the 
flame standoff increases continuously until extinction. 

Figure 7 shows detailed radial profiles of gas temperature (T s ), carbon dioxide mass 
fraction (Tc' 02 ), water vapor mass fraction {Yhio) and fuel vapor reaction rate (wp) 
for the model prediction in Figure 6. This Figure presents three times: immediately 
after ignition (A); midway through the burn. (B); and close to extinction (C). The 
maximum gas temperature decays to approximately 1500 K within 0.1 s after ignition. 
The maximum gas temperature then remains relatively constant until just before 
flame extinction. This behavior is consistent with blowoff extinction where the flame 
extinguishes when the residence time is not sufficient for chemical reaction[15]. While 
the mechanism for extinction is blowoff (as opposed to radiant extinction), radiation loss 
from the flame zone is still important in these flames. The flame temperature and D ext 
for an identical model prediction only neglecting radiant loss are approximately 1600 
K and 0.30 mm, respectively. This emphasizes the potential importance of radiant loss 
even if it is not the dominant mechanism for flame extinction. The model does not 
include broad-band radiative loss (i.e. soot emission), but the flames in the experiments 
are dim and blue so radiative loss from soot is not significant. 

Figure 8 shows the average burning rate constant, k, as a function of ambient 
oxygen mole fraction. The initial droplet size was 1.8 mm and the ambient pressure 
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Figure 6. Numerical model results (compared to experiment) of a Do = 1.7 mm 
droplet in a X 02 = 0.15, P = 120 mmHg ambient. The flame in the experiment exists 
for slightly longer than in the graph, but was nearly invisible and difficult to measure 
(the SiC fiber indicates a flame was present). 


was 190 mmHg for all of the model results. Further, the ignition parameters (location 
and magnitude) were identical for all of the model runs. The initial droplet sizes in the 
experiments ranged from 1.2 to 2.0 mm, and included results in different facilities (GRC 
2.2 and 5.2 second drop towers and the 10 second Japan Microgravity Center), and in 
some cases different ignition energies. The range of experimental conditions explains 
some of the scatter in the experimental data. 

Figure 9 shows comparisons between the model and experiments for average burning 
rate constant (9a) and extinction droplet diameter (9b) as a function of pressure. The 
X 01 and D 0 are 0.15 and 1.8 mm respectively, for the numerical model data. The 
experimental initial droplet diameters ranged from 1.5 to 2.0 mm, and there were 
variations in the facility and ignition conditions. The results for the average burning rate 
constant show reasonable agreement between the model and experiment considering the 
uncertainty in the experimental data and the simplifications in the numerical model. 

The model predictions and experimental data for D ext , however, do not agree 
quantitatively. While the model predicts the experimental D ext at 120 mmHg, we 
note that the pre-exponential factor (.4.) has a value to ensure agreement at this 
particular condition. Because of the influence of the support fiber, the D ext is difficult 
(or impossible) to measure experimentally at pressures higher than 120 mmHg. At 
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Figure 7. Detailed radial profiles of gas temperature ( T g ), carbon dioxide and water 
vapor mass fractions {Yco 2 and Y H 20 , respectively) and fuel vapor reaction rate (wp) 
for three times for the test in Figure 6. 


ambient pressures less than 1207nmiJ<7, the model predicts a much stronger dependence 
of extinction droplet diameter on pressure than that observed experimentally. In fact, 
at pressures below 110 mm Hg, the model shows flame extinction almost immediately 
after ignition (or the model shows no ignition), where the experiments show a spherical 
flame burning for a short time followed by extinction. The discrepancy between the 
model and the experiment in the D ext predictions is due to the simplified single-step 
chemistry in the model. 

Figure 10 shows the predicted influence of ignition on the burning history for three 
tests where the ignition time varied from 1.0 to 1.4 s. Because the droplet history 
for the three tests was nearly identical, the graph displays only one. The nearly 
identical droplet histories means that the burning rate constant is not significantly 
influenced by the ignition time ( k = 0.56, 0.58 and 0.58 mm 2 /s for t ig = 1.0, 1.2 and 
1.4s, respectively). Further, the extinction droplet diameters for the three tests axe 
nearly identical ( D ext = 0.67,0.66 and 0.66 mm for t ig = 1.0, 1.2 and 1.4s, respectively). 
This near independence of k (and k) and D ext on ignition time is in reasonable agreement 
with the small changes in the experiments (Fig. 5). Figure 10 shows that the ignition 
time strongly influences the initial flame size. As the ignition time increases, the initial 
flame position moves further from the droplet. This agrees both qualitatively and 
quantitatively with the experimental results in Figure 5. Later in the droplet lifetime, 
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Figure 8. Comparison of experimental results and numerical model predictions for k 
as a function of X 02 for P = 120 mmHg. 


the model predicts that the flame histories merge and the flame seems to ’forget’ the 
initial condition. The experiments show, however, that the ignition transients persist 
throughout the droplet lifetime. Part of the difference is due to the different pressures 
(90 mmHg in the experiment versus 120 mmHg in the model) which means that the 
droplet lifetimes in the experiments were considerably shorter. The experimental flame 
histories for different ignition times may merge for ambient conditions with longer flame 
lifetimes. 

4. Discussion 

The numerical model accurately predicts the temporal behavior of both the droplet 
and flame. This agreement requires suitable values for the average gas-phase thermo- 
physical and chemical kinetic properties. The fact that the droplet history agrees is 
not surprising, because even the simplest formulation (e.g. [1]) will produce accurate 
estimates of the burning rate constant, as long as reasonable properties are used in 
the formulation. Simplified models, however, do not predict accurate flame sizes. The 
current model provides reasonable agreement for flame size as a function of time. 

One unexpected result of the modeling efforts was importance of matching the 
ignition source term in the model with the igniter characteristics in the experiments. 
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Figure 9. Comparison of experimental results and numerical model predictions for 
(a) k and (b) D ext as a function of P in X02 = 0.15 ambients. 

Specifically, early in the model development, there was a large discrepancy between 
the model and experiment. The model either predicted that the droplet burned to 
completion or extinguished quickly after ignition depending on the value of A. This 
inability to predict extinction disappeared and the model predictions improved when 
the igniter location and energy more closely matched the experiment and the model 
accurately predicted the pre-ignition vaporization behavior. 

One reason for this sensitivity is the degree of pre-ignition heating of the liquid 
phase, particularly for the relatively large droplets in the present study. With the 
ignition region closer to the droplet surface, the droplet ignites quickly but at a relatively 
cold temperature. The rate of heat loss to the cold droplet during combustion then 
significantly affects the extinction behavior of the flame. The heat loss to the droplet is 
directly proportional to the thermal conductivity of the liquid fuel. The current model 
assumes that the energy transport in the liquid is purely diffusional. A number of studies 
(e.g. [19, 20]) exist, however, showing that significant convective flows, driven by surface 
tension gradients, exist inside these relatively large liquid droplets. These surface tension 
gradients can be either thermally or solutally driven, although the later is unlikely in 
the present study. It may be possible to use a modified thermal conductivity [21] to more 
closely approximate the actual thermal transport inside the droplet. We note, however, 
that the presence of the support fiber in the experiments has an unknown influence on 
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Figure 10. Numerical model results for the burning history of a Dq = 1.5 mm droplet 
burning in a P = 120 mmHg, X 02 = 0.15 ambient with three different ignition energies. 
The variation in ignition energy results, similar to the experiments, from changes in 
the ignition time. 


any convective flows inside the droplet. 

The sensitivity to ignition may occur because combustion takes place in ambients 
that are close to the flammability limit. The current model predicts a limiting Y 02 of 
approximately 0.14 that can support combustion for the droplet size and pressure range 
in this study. The limit predicted by the model is nearly identical to the experimental 
limit observed by Easton[22], 

The transient behavior of the flame standoff in the present experiments suggests 
the potential importance of fuel vapor accumulation [24, 25] between the droplet and 
the flame. When the fuel vaporization rate is greater than the total fuel consumption 
rate at the flame, fuel can accumulate between the droplet and the flame. Law and 
co-workers [24] showed that this phenomena can explain the experimentally observed 
growth in the flame standoff ratio with time (especially early in the droplet lifetime). 
The ratio (5) of the fuel vapor consumption rate to the fuel vaporization rate is easy to 
evaluate in the present model. 

2 D 2 r°° 

6 = — — / f 2 wpdr ( 10 ) 

pL K J 1 

The wp in Eq. 12 is the radial variation of the fuel vapor consumption rate such as that 
in Fig. 7. 
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He’s [25] analysis indicates that the motion of the flame front is important and 
the fuel vapor consumption rate at the flame can differ from the fuel vaporization rate 
for the ambient conditions in the present study. In this case, their model predicts a 
flame size that increases initially with time, reaches a maximum and then decreases in 
response to the shrinking droplet, (^he^nodel further predicts that the flame standoff 
ratio will increase continuously with time. This is in qualitative agreement with the 
flame behavior in the present experiments. 


Jhju 



Figure 11. Fuel vapor consumption to fuel vaporization rate ratio (<5) as a function 
of non-dimensional time for different X 02 ambients (P = IQOmmHg). The flame data 
are for the X 02 = 0.18 ambient. 


Figure 11 shows the model prediction of the temporal variation in <5 with time for 
different X 02 ambients ( P = 190 mmHg). For X 02 greater than 0.16, 6 is greater than 
0.9 for nearly all of the droplet lifetime, and gradually increases throughout the burn 
approaching nearly 1.0 before decreasing abruptly near fuel droplet burnout. Figure 11 
also shows the temporal variation of D f and Dj/D for the Y 0 2 = 0.18 case. The results 
clearly show that 6 is close to 1.0 even though Df/D increases continuously until very 
close to droplet burnout. 

Figure 11 further shows that when X 02 = 0.16, 5 increases quickly to 0.8, then 
gradually increases to 0.9 before decreasing (more gradually than the higher X 02 cases) 
near fuel droplet burnout. The temporal behavior of 5 is very different in the X 02 = 0.15 
ambient, with the 5 less than 0.6 for almost the entire droplet lifetime. A detailed 
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examination of the w F profiles, however, shows this decrease is, at least in part, due to 
fuel leakage through the flame front. 

5. Conclusions 

This paper presents experiments involving the combustion of single decane droplets 
in microgravity. The experiments are over a range of ambient pressures and oxygen 
mole fractions, initial droplet sizes and ignition energies. The data are useful for the 
validation of detailed numerical models of droplet combustion which in turn can form 
the basis of sub-models for spray combustion models. The data regarding extinction 
droplet sizes enable the evaluation and validation of kinetic mechanisms in a simplified, 
computationally tractable geometry. 

The numerical model accurately predicts the burning behavior of the droplet, both 
the droplet and flame history. This good agreement requires accurate modeling of 
the ignition which in turn requires accurate modeling of the pre-ignition vaporization 
behavior. Both the experiments and model show transient flame behavior which suggests 
the importance of fuel vapor accumulation between the droplet and flame. The model 
results show, however, that the ratio of the mass rate of fuel consumption (at the flame) 
to fuel vaporization is nearly 1.0 suggesting that fuel vapor accumulation alone cannot 
explain the observed transient flame behavior. 
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It-* 


Nomenclature 


A 

A i P 

A 

C p 

Vi 

D 

D 0 

D f 

Da 


mean absorption coefficient 

Planck mean absorption coefficient for species i 

Pre-exponential factor (8 x 10 14 cm 3 /g/s ) 

gas-phase specific heat (0.362 cal/g/K) 

binary diffusion coefficient of species i 

droplet diameter 

initial droplet diameter 

flame diameter 

Damkohler number 
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E 

k 

k 

C 

Lei 

Mi 

M 

Pi 

Pv 

Q 

Qr 

Qrs 

Qi 

r 

f 

r 

t 

R 

n 

T 

T b 

u 

Wp 

X 

Y 

a 

oil 

A 

Al 

$ 

<t> 

p 

PL 

Vi 


activation energy (3 x 10 4 cal /mole) 
instantaneous burning rate constant 
average burning rate constant 

latent enthalpy of vaporization (varies with T, 70 cal/g at T b ) 

Lewis number of species i (2, 1.1, 1.4 and 0.83 for F, 0 2> C0 2 and H 2 0, respectively) 
molar mass of species i 
average molar mass 
partial pressure of species i 

equilibrium fuel vapor pressure at surface temperature 
heat release per unit mass of fuel (11472 cal/g) 
gas phase radiative energy loss (cal/ s/ cm 3 ) 
surface radiative energy loss (cal/ s/ cm 2 ) 
ignition energy (cal/ s/ cm 3 ) 
radial coordinate 

non-dimensional radial coordinate r — % 
transformed radial coordinate f — tan f 
time 

radius of the droplet 
universal gas constant 
temperature 

boiling temperature of the fuel 
gas velocity 

fuel vapor reaction rate per unit volume (g/ s/cm 3 ) 
mole fraction 
mass fraction 

gas-phase thermal diffusivity 
liquid-phase thermal diffusivity 

gas-phase thermal conductivity (1.05 x 10 -4 cal/K/s/cm) 
liquid-phase thermal conductivity (2.82 x 10 -4 cal/K/s/cm) 

^ PocCp Too 

velocity potential, u T = 
gas density 
liquid fuel density 

stoichiometric coefficient for species i 
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